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INTRODUCTION 


raining information about cirrus cloud microphysics requires development of remote sensing techniques. 
In an earlier paper, Matrosov et al. (1992) proposed a method to estimate ice water path (IWP)(ne„ vertically 
integrated ice mass content IMQ and characteristic particle size averaged through the cloud from combined ground- 
based measurements of radar reflectivities and IR brightness temperatures of the downwelling diermal radiation in 
the transparency region of 10-12 pm. For some applications, the vertically averaged characteristic particle sizes and 
IWP could be the appropriate information to use. However, vertical profiles of cloud microphysical parameters can 

provide a better understanding of cloud structure and development. . 

Here we describe a further development of the previous method by Matrosov et al. (1992) for retrieving 
vertical profiles of cirrus particle sizes and IMC rather than their vertically averaged values. In addition to 
measurements of radar reflectivities, the measurements of Doppler velocities are used m the new method. This 
provides us with two vertical profiles of measurements to infer two vertical profiles of unknowns, i.e., partic e 
characteristic sizes and IMC. Simultaneous measurements of the DR. brightness temperatures are still needed to reso ve 
an ambiguity in particle size - fall velocity relationships. 


2. THEORETICAL CONSIDERATIONS 

Doppler velocities, V m . measured with the vertically pointed radar antenna represent the sum of reflectivity- 
weighted particle fall velocities, V f , and of vertical air motion, V a : 

( 1 ) 


v m = v,+ v a . 

To extract values of V. from measured quantities V m , we used an approach based on the assumption that, over time 
intervals of 1 - 2 hours, the mean vertical air motion in a cloud is small compared to fall velocities of ice particles. 
Thus, the observed average value of Doppler velocity closely approximates F/Orr and Kropfli 1993). Averages were 
performed for 1 dB reflectivity intervals to find reflectivity- fall velocity relationships at different range gates to get 
estimates of V, at the same time intervals as those of reflectivity. Our experience with Doppler velocity data 
processing shows that changing the averaging time from 1 hour to 2 hours usually does not cause significant 
variations in estimated values of fall velocities which demonstrates the soundness of this approach. 

The largest dimensions of ice particles in cirrus clouds usually do not exceed - 2 mm (Dowling and 
Radke 1990), which is still within the Rayleigh scattering regime for radar frequencies up to the Ka-band 1 (Yeh et 
al 1982) In this regime, backscattering cross sections of individual particles increase with the sixth power of particle 
characteristic size. After integrating these cross sections with the particle size distribution, we get. for the radar 
reflectivity with respect to ice, Z,, 


Z,=f,CD m 


( 2 ) 
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where C and D m are particle concentration and median diameters, respectively, and the coefficient /, depends on the 
particle size distribution. Hereafter, we assume that this distribution is the gamma function of the first order. This 
function adequately describes many experimental size spectra of cirrus cloud particles (Kosarev and Mazin 1989). 
For the first order gamma function, /, * 0.486 if Z i is in cm^. C is in cm' 3 , and D m is in cm. Equation (1) implies 
that particles are spherical. Effects of nonsphericity, however, can be approximately taken into account by tuning 
the coefficient /,. 

Integrating over the distribution, we can obtain also the expression for EMC as follows: 

IMC = f 2 CD m \ (3) 

where f 2 ~ 0.111 if the size distribution is the first order gamma function and EMC is in g cm 3 . 

From (2) and (3) one can see that finding vertical profiles of particle sizes and IMC is equivalent to having 
estimations of C and D m at each radar range gate. Equation (2) provides one vertically resolved relationship between 
the measured values (i.e. Z) and the unknowns. A second relationship should connect these unknowns and the 
particle fall velocity estimates. 

Particle fall velocity V f is the reflectivity-weighted velocity of individual particles in the radar resolution 
volume, v f : 
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V f = / v J N(D)D 6 dDIZ i , 
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(4) 


where N(D) is the size distribution function. 

Experimental studies of the fall velocities of individual ice crystals demonstrate that fall velocity-size 
relationships can be approximated by the power law function (Pruppacher and Klett 1978): 

ss AD 8 , (5) 

where A and B are constant for a particular crystal shape. According to the data presented by Pruppacher and Klett 
(1978), the coefficient B generally varies from 0.8 to 1.3. The coefficient A shows much greater variations, up to 
two orders of magnitude. 

Integrating according to (4) gives the following expression for V f estimates: 

V f = Af 3 D m B . (6) 

The coefficient f 3 depends slightly on D m , but this dependence is negligibly small. The value of f s depends also on 
B . For the first-order gamma function and B=Uf 3 » 1.71 if V f is in cm s' 1 . 

Equation (6) provides the second vertically resolved relationship between the measurables and unknowns. 
Given relatively low variations of 5, we can reasonably assume that B * 1. However, possible large variations of 
A indicate that the value of this coefficient has to be estimated from at least one additional measurement. 

A measurement for estimating A is obtained from the brightness temperature of cloud downwelling thermal 
radiation. Matrosov et al. (1992) proposed a technique to infer cloud optical thickness due to extinction from 
brightness temperature measurements in the ER "window” region by a narrow-band ER radiometer. In this region, 
cirrus particle size factors ( kD/X ) are large compared to the wavelength, and we can assume that particle extinction 
efficiency is close to 2. In this case, the cloud extinction coefficient a and the optical thickness t can be expressed 
as follows: 


(7) 

( 8 ) 

where the summation is with respect to range gates A hj and f 4 « 0.432 if a is in cm' 1 . 


and 


a =f£D m \ 

T = E VJCDJM. 
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At ^ firs, step of 0* retrieval cloud tnoophyskal parameter fa »J « «** some **£ 

value for A. Then, from estimates of fall velocities V, we retrieve values of P^'e ' 3 f < » 8 

gam. Knowing D. allows calculations of panicle ^cenuanons T ftom m radiomeBr 

value of optical thickness t is calculated using (8) and compare measured values of optical 

measurement for the considered radar beam. From die ratio of “of t would cohScide. 

thickness we find a corrected value of A for which the newly ca c . . ■ concentrations, and IMC 

The corrected value of A is then used to calculate final vertical profiles of particles sizes, concentrauo 

using (6). (1), and (3). 

3. experimental EXAMPLE 

We use experimental data obtained during Ore FIRE-n experiment to Ulustralt i the 
of the priority dams was November 26. when a slowly with 

S '• K ST W« " S^cmwav e 

a wavelength band from 9.9b to ii_ 3 iun on f, GHz') to account for the atmospheric transparency and 

therm^^ad^rtion^vhen reriev^ optical thickness values from IR brightness temperatures usmg the technique 

is the following: r e = 0.46 D m . 



Equivalent Median Diameters of Cirrus Particles 26-NOV-1991 


D^m) 


400 0 
3500 
300 0 
250 0 
200.0 
150 0 
100 0 
50 00 



Figure 1. Time height cross 


section of cirrus particle median diameters observed on November 26, 1991 
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Figure 2 shows a time height cross section for retrieved values of particle concentrations. The highest 
particle number densities are in the upper part of the cloud. However, the sizes of these particles are small. As 



Figure 2. Time height cross section of cimis particle concentrations observed on November 26, 1991. 

The size information obtained using the radar/IR radiometer technique is in general agreement with the 
lidar/radar te chni que and data of direct measurements (Intrieri et al. 1993). Knowing characteristic sizes, 
concentrations and te rmin al velocities of particles one can easily calculate some other parameters of cloud 
microphysics such as vertical profiles of ice mass content and ice mass flux. Integrating vertical profiles of ice mass 
content provides values of ice water path which is important for longwave cloud feedback (Ebert and Curry. 1992). 
Our future plans include retrieving and studying information about aforementioned parameters and estimating possible 
retrieval errors. 
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